Measurement of Ultrafast Carrier Dynamics in Epitaxial Graphene 
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Using ultrafast optical pump-probe spectroscopy, we have measured carrier relaxation times in 
epitaxial graphene layers grown on SiC wafers. We find two distinct time scales associated with the 
relaxation of nonequilibrium photogenerated carriers. An initial fast relaxation transient in the 70- 
120 fs range is followed by a slower relaxation process in the 0.4-1.7 ps range. The slower relaxation 
time is found to be inversely proportional to the degree of crystalline disorder in the graphene 
layers as measured by Raman spectroscopy. We relate the measured fast and slow time constants 
to carrier-carrier and carrier-phonon intraband and interband scattering processes in graphene. 



Graphene is a single two dimensional (2D) atomic layer 
of carbon atoms forming a dense honeycomb crystal lat- 
ticed. It is a zero-bandgap semiconductor with a linear 
energy dispersion relation for both electrons and holes^. 
The unusual electronic and optical properties of graphene 
have generated interest in both experimental and theo- 
retical arena a 2 ' 3 ' 4 ' 5 ' 6 . The high mobility of electrons in 
graphene has prompted a large number of investigations 
into graphene based high speed electronic devices, such 
as field-effect transistors, pn-junction diodes and tran- 
sistors, and terahertz oscillators, and also into low noise 
electronic sensor s 4 ^^ ' 1 1 1 1 . 

The simplest way of obtaining graphene layers is via 
micromechanical cleaving (exfoliation) of bulk graphite 
followed by careful selection of monolayers by using op- 
tical, atomic force, or scanning electron microscopes^. 
Although this technique results in relatively high quality 
films, it might not be suitable for large scale production. 
Recently, epitaxial growth of graphene by thermal de- 
composition of SiC surface at high temperatures has been 
investigated as a promising alternative for large scale pro- 
duction 6,12 . This technique can provide anywhere from 
a few monolayers of graphene to several (> 50) layers 
on the surface of a SiC wafer. Graphene layers grown 
by this technique have demonstrated structural and elec- 
tronic properties similar to those of graphene layers ob- 
tained by micromechanical cleaving techniques, includ- 
ing the massless Dirac-like energy dispersion relation for 
electrons and holes and carrier mobilities in the few tens 
of thousand cm 2 /V-s range 6 - ! 13 ' 14 . In addition, the elec- 
tronic as well phononic properties of epitaxially grown 
graphene multilayers have been found to be different from 
those of bulk graphite and similar to those of a graphene 
monolayer indicating that the electrons and phonons in 
different layers in epitaxially grown graphene are uncou- 
ple d 16 ' 17 . This observed difference in the properties of 
epitaxial graphene and bulk graphite has been attributed 
to a different stacking scheme for carbon atom layers in 
epitaxial graphene compared to bulk graphite 1 ^-. Epi- 
taxial growth of graphene on SiC provides a technique to 
obtain large area multilayers that can be used for studies, 
such as ultrafast optical spectroscopy, that are difficult 
to conduct on monolayers. 

In this paper, we present results from measurements 
of the ultrafast dynamics of photoexcited carriers in 



graphene for the first time. Ultrafast studies of car- 
rier dynamics in other forms of carbon, such as car- 
bon nanotubes and bulk graphite, have been reported 
in the pas t 18 ' 19 ' 20 . The results presented in this paper 
are relevant for understanding carrier intraband and in- 
terband scattering mechanisms, and the corresponding 
rates, in graphene and their impact on proposed and 
demonstrated graphene based electronic and optical de- 
vices*>LS&i&. 

The graphene samples used in this work were all epi- 
taxially grown on the carbon face of semi-insulating 6H- 
SiC wafers using the techniques that have been reported 
in detail previously 4 ^. Samples A, B and C were grown at 
temperatures varying from 1400 °C to 1600 °C and pres- 
sures of 2 — 7 x 10~ 6 torr. Micro-Raman spectroscopy of 
all samples showed a single-resonant G peak close to 1580 
cm" 1 , a double-resonant D' peak close to 2700 cm -1 , and 
also a relatively low intensity double-resonant D peak 
near 1350 cm" 115 ' 16 . The D peak is not allowed in per- 
fect graphene layers since it requires an elastic scattering 
process, which is made possible by disorder, to satisfy 
momentum conservation 16 . The presence of the D peak 
therefore indicates the presence of disorder in the sam- 
ples. A discussion of the observed intensity of the D 
peak and its correlation with the measured carrier re- 
laxation times is presented later in this paper. Fourier 
Transform Infrared (FTIR) spectroscopy of the samples 
revealed a fiat absorption profile in the entire 2.5-25 /im 
range, which is consistent with the massless Dirac-like 
energy dispersion of electrons and holes (see Eqs.[T]and 
[2] below). The number of graphene layers in samples A, 
B, and C were estimated from FTIR and XPS spectro- 
scopies (using the Thickogram method 24 -) to be 6, 12 and 
37, respectively, with less than 5% error. 

A Tksapphire mode-locked laser with 86 MHz pulse 
repetition rate, 780 nm center wavelength, and ~85 fs 
pulse width was used for time-resolved pump-probe spec- 
troscopy of the graphene samples. Pump pulses with 
energies between 3-15 nJ were used to generate photoex- 
cited carriers, while weak probe pulses with energies be- 
tween 30-100 pJ were used to measure the changes in 
the transmittivity of the samples at various delays of the 
probe pulses with respect to the pump pulses. The an- 
gle of incidence of the pump and probe beams were 0° 
and 15° respectively. The pump and the probe were fo- 



2 




t = 0.1 2 ps ( a ) 
t = 1 .4 ps Sample A 

2 Pump pulse energy = 12.5 nj 



(b) 



Sample B 



Pump pulse energy = 13.4 nj 



(c) 

Sample B 
Pump pulse energy = 8.8 nj 



(d) 

_ 0.4 ps Sample C 

2 Pump pulse energy = 9.9 nj 



2 3 
Delay (ps) 



FIG. 1: Measured transmittivity transients for (a) sample 
A (b) sample B (c) sample B with different pump power, 
and (d) sample C. The dark solid lines with markers are the 
experimental data and the light solid lines without markers 
are analytical fits to the data using exponentials with time 
constants ti and T2. The transients in (b) and (c) show that 
the slower time constant T2 does not change much as the pump 
energy is varied. 



cused to a spot size of about 100 /im. The polarization 
of the probe was rotated by 90 degrees with respect to 
the pump and a polarizer was used to eliminate scattered 
pump light going in the direction of the probe beam. The 
probe beam was passed through a 50 ^m spatial filter for 
further removal of the scattered pump light. The pump 
and probe beams were both modulated at two different 
frequencies near ^3 KHz, and changes in the intensity 
of the probe pulses at the sum of these two frequencies 
were measured with a lockin amplifier. 

Photon interband absorption in graphene at optical 
and near-infrared frequencies is given by the optical con- 
ductivity a(ui)r^ : 



cr(w) 



(1) 



where, f c (E) and f v (E) are the probabilities for the oc- 
cupation of an energy level with energy E in the conduc- 
tion and valence bands, respectively. The only frequency 
dependence of cr(w) comes from the carrier distribution 
functions. Matching the optical boundary conditions at 
the air/graphene/SiC interfaces, the optical transmission 
t{ui) through N graphene layers on a SiC wafer (normal- 
ized to the transmission through a plain SiC wafer) can 




FIG. 2: (a)Band structure of graphene showing an intrin- 
sic population of electrons and holes near the Dirac point. 
Optical excitation is indicated by the arrow. (b)The non- 
equilibrium distribution of photoexcited carriers account for 
the initial rise in transmittivity. (c) The carriers equilibrate 
among themselves through carrier-carrier scattering on a time 
scale given by n resulting in a hot carrier distribution, (d) 
Subsequent cooling and decay of the hot distribution through 
carrier-phonon scattering (and possibly electron-hole recom- 
bination) occurs on a time scale given by T2. 



be written as, 
t(w) 



1 + Na{u))y/fj, /e /(l + nsic) 



(2) 



where, nsic is the refractive index of SiC. The above ex- 
pression can be used to estimate that pump pulses with 
energies in the range indicated above generate electron 
and hole densities in the 3 x 10 11 -10 12 cm' 2 range. The 
photogenerated carrier densities are larger than the in- 
trinsic electron and hole densities of ~ 8 x 10 10 cm~ 2 in 
graphene at room temperature. 

Figure 1 shows the measured transmittivity transients 
for different graphene samples. The figure shows the time 
dependent change AT in the transmittivity normalized 
to the transmittivity in the absence of the pump pulse. 
Transmittivity increases sharply just after photoexcita- 
tion. The recovery of the transmittivity exhibits two dis- 
tinct time scales; an initial fast relaxation time n in the 
70-120 fs range followed by a slower relaxation time 
in the 0.4-1.7 ps range. These time constants have been 
extracted by analytical fits to the data using decaying ex- 
ponentials. It should be noted here that the faster time 
T\ is of the order of the pulse width and is therefore not 
accurately resolvable. 

A simple model incorporating band-filling effects to- 
gether with intraband carrier-carrier and carrier-phonon 
scattering can be used to explain the observed trans- 
mittivity transients. Figure 2 is a schematic represen- 
tation of this model. As shown in Eq. [TJ optical inter- 
band absorption in graphene is proportional to the dif- 
ference between the occupancies of the conduction and 
valence bands at energies equal to Tiw/2 (measured from 
the Dirac point). Photogeneration of carriers by the 
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pump pulse reduces this difference and causes the ini- 
tial increase in the transmittivity observed in Figure 1. 
Immediately after photoexcitation, the non-equilibrium 
carrier distribution broadens and also equilibrates with 
the intrinsic carrier population through carrier-carrier 
scattering. This process results in the initial fast relax- 
ation of the transmittivity. The observed fast relaxation 
times (rx ) are consistent with the theoretically predicted 
carrier-carrier intraband scattering rates in graphene by 
S. Das Sarma et. ali^S. As a result of carrier-carrier 
scattering the photogenerated carriers are expected to 
equilibrate among themselves and reach a Fermi-Dirac- 
like distribution with a temperature much higher than 
the lattice temperature. The observed slower time con- 
stant (j2) of the transmittivity decay could be attributed 
to the subsequent thermalization of the carriers with the 
lattice through carrier-phonon intraband scattering. 
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FIG. 3: The slower relaxation time T2 is plotted versus the 
ratio of the intensity of the Raman G and D peaks for samples 
A, B and C. This ratio is a measure of the crystal coherence 
length. Larger crystal disorder (smaller coherence length) re- 
sults in shorter relaxation times. 

Electron-hole recombination processes could also con- 
tribute to the slow decay of the transmittivity. The 
dominant mechanisms for electron-hole recombination in 
graphene are not yet well understood. Electron-hole re- 
combination due to Auger scattering in graphene was 
analyzed and carrier density dependent lifetimes of the 
order of a few picoseconds for electron-hole densities in 
the I0 11 -I0 12 cm -2 range were predicted by F. Rana 24 . 
However, in our experiments varying the pump pulse en- 
ergies in the 3 nJ to 15 nJ range to vary the photo- 
generated carrier densities did not lead to any signifi- 
cant changes in the measured values of t%. For example, 
Figures 1(b) and 1(c) show the measured transmittiv- 
ity transients for two different pump pulse energies for 
sample B. These results indicate that the dominant con- 
tribution to T2 comes from a process that is independent 
of the carrier density and is likely carrier-phonon rather 
than carrier-carrier scattering. However, electron-hole 
recombination due to carrier-phonon interband scatter- 
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FIG. 4: Maximum increase of the transmittivity after pho- 
toexcitation is plotted for various pump energies for sample B. 
Linear relation between the maximum transmittivity change 
and the the pump pulse energy agrees with the linear ab- 
sorption and band filling model and rules out nonlinear two- 
photon absorption. 



ing cannot be ruled out. It needs to be pointed out here 
that carrier-phonon scattering cannot also be ruled out 
as a contributor to the fast relaxation time n from our 
measurement results. Pump pulse energy dependence, 
and therefore carrier density dependence, of the time t\ 
could not be reliably extracted from the measurements 
since, as already mentioned above, the observed values of 
Ti were close to the pulse width used in the experiments. 

It has been shown that carrier-phonon deformation po- 
tential scattering rates in semiconductor nanostructures 
are enhanced in the presence of disorder—. In graphene, 
the intensity of the double-resonant D peak (near 1350 
cm -1 ) in the raman spectrum can be used as a measure 
of crystalline disorder since this peak is absent in per- 
fect graphene layer o 16 i 26 . The ratio of the intensities, Iq 
and Id, of the G and D peaks, respectively, in the ra- 
man spectrum has been shown to be proportional to the 
crystal coherence length^. Thus, one could expect the 
measured time constant to scale with the ratio Ig/Id- 
Figure 3 shows the measured values of T2 plotted vs the 
measured values of the ratio Ig/Id for the three samples. 
Figure 3 shows that T2 oc Ig/Id and therefore T2 is pro- 
portional to the coherence length of the crystal. These 
results also support electron-phonon scattering as being 
the dominant mechanism contributing to t%. 

Figure 4 shows the peak value (normalized) of the mea- 
sured transmittivity change AT/T plotted as a function 
of the pump pulse energy for sample B. As expected from 
linear absorption and final state filling arguments, the 
maximum value of AT/T is proportional to the pump 
pulse energy. This data also rules out any significant role 
played by nonlinear two-photon absorption in the trans- 
mittivity transients. Complete saturation or bleaching 
of the absorption is not observed for the range of pump 
pulse energies used in the experiments. From the den- 
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sity of states of graphene it follows that the maximum 
electron density in an energy interval AE centered at E 
is 2E 'AE I '-nhf 'v 2 , where v is the velocity of carriers in 
graphene and equals 10 s cm/A Putting E = Tilo/2 and 
AE = TiAuj/2, where Acu is the optical bandwidth of the 
pump pulse (~10 nm), the maximum electron (or hole) 
density comes out to be 1.1 x 10 12 cm~ 2 for 780 nm pump 
center wavelength. The fact that no bleaching of the ab- 
sorption is observed even for pulse energies large enough 
to excite electron (and hole) densities close to 10 12 cm~ 2 
is because of the fast relaxation time t\ that is of the 
order of the pump pulse width. Graphene therefore has 
potential for use as a fast saturable absorber for gener- 
ating high energy short pulses from modelocked lasers^. 

In conclusion, we have measured ultrafast carrier re- 
laxation rates in epitaxially grown graphene layers on 
SiC. We observe two distinct time scales associated with 
the relaxation dynamics of photogenerated carriers. The 
observed time scales are comparable to those observed in 
other related forms of carbons, such as highly ordered py- 
rolytic graphite (HOPG) and single-walled carbon nan- 
otube a 18 ' 19 ? 20 . Our measurements indicate the separate 
roles played by carrier-carrier and carrier-phonon scat- 
tering in relaxing photogenerated carriers. 

More work is needed to investigate the role of carrier- 



phonon scattering in the fast relaxation time T\ , the role 
of electron-hole recombination in the slow relaxation time 
T2, and the value of T2 in the limit of disorder- free epitax- 
ial graphene layers. It has been recently pointed out that 
the first few carbon layers in epitaxially grown graphene 
acquire a bandgap as a result of interaction with the 
atoms in the SiC substrate that breaks the symmetry 
between the A and B atoms in the graphene lattice^. 
The effect of bandgap on ultrafast intraband and in- 
terband carrier dynamics is not clear. Also, in multi- 
layer graphene structures the optical response is likely 
to be dominated by the large number of layers that are 
not close to the substrate and do not have a bandgap. 
Pump-probe experiments with fewer monolayers of epi- 
taxial graphene would be needed to explore the effects of 
bandgap on carrier dynamics. 
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